Displacements of 180° domain walls in electroded ferroelectric 

single crystals: 
the effect of surface layers on restoring force 

A. Kopal/'Q P. Mokry/'l J. Fousek, 1 * 2 and T. Bahmk 1 

1 Dept. of Physics, International Center for Piezoelectric Research, 
Technical University Liberec, Liberec 1, 461 17 Czech Republic 
2 Materials Research Laboratory, The Pennsylvania 
State University, State College, PA 16801, USA 
(Dated: February 6, 2008) 

Abstract 

Macroscopic properties of ferroelectric samples, including those in form of thin films, are, to 
large extent, influenced by their domain structure. In this paper the free energy is calculated for 
a plate-like sample composed of nonferroelectric surface layers and ferroelectric central part with 
antiparallel domains. The sample is provided with electrodes with a defined potential difference. 
The effect of applied field and its small changes on the resulting domain structure is discussed. 
This makes it possible to determine the restoring force acting on domain walls which codetermines 
dielectric and piezoelectric properties of the sample. Calculations of the potential and free energy 
take into account interactions of opposite surfaces and are applicable also to thin films. 
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I. INTRODUCTION 



It is known that samples of ferroelectric single crystals often posses a surface layer whose 
properties differ from those of the bulk. It may be a layer produced during the growth 
of a crystalline plate or produced during the preparation of a plate-like sample. Many 
observations gave evidence to the fact that such a layer is either nonferroelectric or does 
not take part in the switching process of the internal part; in any case its permittivity is 
believed to differ from that of a homogeneous sample in the ferroelectric phase. Its existence 
,s expected to great* influence _p,c ptopett.es of bulk satnpl JflQQ as we„ as 
of thin films. [51] In this paper two such consequences are investigated. First we reconsider 
the problem of equilibrium domain structure in a ferroelectric sample possessing a surface 
layer, previously discussed by Bjorkstam and Oettel.|6i] Second, we evaluate the restoring 
force acting on 180° domain walls due to the layer; this will make it possible to estimate the 
extrinsic contributions to permittivity, piezoelectric coefficients and elastic compliances of a 
ferroelectric sample. Investigations of crystals of the KDP family revealed the existence of 

un 

huge wall contributions to these properties. [3, |8J] 

In previous papers on a related subject j^j], depolarizing field was considered as the source 
of energy which slows down the motion of a single domain wall in a dc electric field, as the 
wall departs from its original position by substantial distances. In contrast to such models 
we investigate very small deviations of walls forming a regular domain pattern. 



II. GEOMETRY, VARIABLES AND ENERGY OF THE SYSTEM 

We consider a plate-like electroded sample of infinite area with major surfaces perpendic- 
ular to the ferroelectric axis z. Central ferroelectric part with antiparallel domains (2.) is 
separated from the electrodes (0.), (4.) by nonferroelectric layers (1.), (3.) (see Fig.[T]). The 
spatial distribution of the electric field E is determined by the applied potential difference 
and by the bound charge divPo on the boundary of ferroelectric mate- 
rial, where Po stands for spontaneous polarization. Geometrical, electrical and material 
parameters of the system are shown in Fig.[TJ 

We further introduce the symbols 

c = ^/e a /e c , g = y/e a e c , 
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FIG. 1: Geometry of the model 
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The ferroelectric material itself is approximated by the equation of state 
D z = E$E C E Z + P , 

where Pq is the spontaneous polarization along the ferroelectric axis. This linear approxi- 
mation limits the validity of our calculations to the temperature region not very close below 
the transition temperature T c . Domain walls are assumed to have surface energy density a w 
and zero thickness. 

The Gibbs electric energy of the system includes the domain wall energy, the electrostatic 
energy whose density is (1/2) JS • (D — Pq) and the work performed by external electric 
sources —VQ, where Q is the charge on positive electrode. First, Laplace equations have to 
be solved for electric potentials in the bulk and in the surface layers, fulfilling the requirement 
of potential continuity as well as conditions of continuity of normal components of D and 
tangential components of E. A rather cumbersome calculation leads to the following formula 
for Gibbs electric energy per unit area of the system (in Jm -2 ): 



4P|t 



E 



sin 2 (rwr(l + A)/2) 



G = -a R >p A P « ABt - 2 toe d V(l + B) 

it w 2e (l + B)(e d + Be c ) e TT 2 R(l + B) ^ n 3 (coth nBR + #coth nRc) 
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The first term represents domain wall contribution while the last one is the depolarization 
energy. In the second term we recognize the effect of layers (1.) and (3.) and of the applied 
voltage. 

Let us compare expression with formulae deduced and used in previous papers. For V = 
and s — > oo, Ed = 1 and A = the system goes over into an isolated ferroelectric plate with 
"neutral" domain structure, placed in vacuum. In this case the equation (pQ) reduces to 
the expression given by Kopal et al. 9] for ferroelectric plates of finite thickness in which 
the interaction of the two surfaces is accounted for. If the plate is thick this interaction 
can be neglected and Eq. ([T|) simplifies to the classical formula of Mitsui and Furuichi 10J 
(cf. Eq. (9) in Ref. [9j) which is often used to determine the value a w from the observed 
width of domain patterns. Finally, in the limit of V = our formula (0Q) should converge 
to the expression deduced by Bjorkstam and Oettel. [6J In fact this is not the case and it 
appears that the electric displacement as expressed in Ref. Q| does not satisfy all boundary 
conditions. 



III. EQUILIBRIUM DOMAIN STRUCTURE FOR V = 



If the system is short-circuited, the equilibrium domain pattern is symmetric, i.e. A eq = 0. 
The shape factor R® q and from it also the value of a° q can be found by numerical methods. 
As an example, the full lines in Fig. [2] shows the a eq (P) dependence at constant hid and 
different values of potential difference V. As it was shown in the previous paper 9], the 
critical thickness h crit can be defined by the formula 

^-crit — 4 £o<y w y/e a e c / Pq 

so that for h ^> h CTit the interaction energy of sample surfaces can be neglected. Then the 
minimum energy occurs for 

1/2 

(2) 



eq 



1 



7T 



8S (1 + B) 



y/t. 



and the dependence ctg q (£) in this approximation is shown by a dashed line with V = V 
in the Fig.0 The approximate results in the Fig. [2] with V = 250 V and V = 500 V are 
based on Eqs. and (jl]). In these numerical calculations we have used the following values 
which roughly apply to crystals of RbH 2 P04 below the transition temperature: Pq = 5.7 x 




FIG. 2: Exact numerical (full lines) and approximate (dashed 
lines) results for a eq (t) at different values of potential difference V 
and B = 0.02. 



1(T 2 C m~ 2 , e a = 10, e c = 100, e d = 10, h ait = 5.4 x 10~ 8 m. The value a w = 5 x lO" 3 J m" 2 
is often considered typical for ferroelectrics. 



IV. RESPONSE OF DOMAIN STRUCTUTE TO EXTERNAL ELECTRIC FIELD 

When an external potential difference V is applied, the asymmetry parameter becomes 
nonzero and at the same time the period a = (a + + a_)/2 changes. Both these quantities can 
be found again by determining the minimum of G given by Eq. (pQ) numerically. Full lines 
in Fig. [3] demonstrate both these dependencies for the following numerical values: B = 0.02, 
t = 10 4 t cr it = 10 4 (1 + B) h CTit . Starting from a certain applied voltage the domain spacing 
a grows very fast with increasing V. To discuss the macroscopic properties of the sample 
such as permittivity, the dependence A eq (V) is more important. We recognize that in a 
considerable region of the applied voltage this dependence is almost linear. 

The A eq (V), a eq (V) resp. R eq (V) dependence can be approximated by an explicit formula 
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FIG. 3: Exact numerical (full lines) and approximate (dashed 



lines) results for a e(l (V) and A eq (V) at B = 0.02 and t = 10 4 t. 



crit- 



if the following inequalities are satisfied: A <C 1, BR ^> 1 and Rc ^> 1. Then it holds 



A eq {V)=e e d (l + B) 



V/t 



B- 



2\og2{e d + Be c ) 



(3) 



R eq (V) = R eq ^j 1 - ^ Al^V) . (4) 

These approximations, shown in Fig.[3]by dashed lines, are based on the limit of the sum in 
Eq.©. 



V. EXTRINSIC PERMITTIVITY 



The nonzero value of the asymmetry parameter means that an extra bound charge is 
deposited on the electrodes due to the domain wall shifts and this in turn represents an 
increase of effective permittivity of the whole system crystal plus both surface layers. The 
increase of the electrostatic energy when domain walls leave their original equilibrium posi- 
tions for E ^ serves as the source of a restoring force when field is again reduced to zero. 
The calculations show that the effective permittivity e e ff defined by the total capacitance 
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per unit area C = e e e s/t equals 

r i -i 

R R 9 1n9h J. Rc /c.1 

(5) 



g _ £ C (1 + B) 1 + B ' 

eff 1 + 5£ c /£d 1 + Be c /e d 



g 21n2(l + £e c /e d ) 
e d i?o q (e d + y/£ a e c ) 

In this formula the first term on the right-hand side represents the intrinsic part of permit- 
tivity, given by linear dielectric response of the sample and of the surface layers when domain 
walls are kept at rest. The second term is the contribution of domain walls displacement to 
effective permittivity, often referred to as extrinsic part of permittivity 

VI. DISCUSSION 

Numerous data are available on domain wall contributions to permittivity in single crys- 
tals of ferroelectrics and also on extrinsic contributions to piezoelectric coefficients in ferro- 
electrics which are simultaneously ferroelastic. Our calculations indicate that depolarizing 
energy can be an effective source of restoring force whose existence is a condition for such 
contributions. In fact since the model assumes a regular system of planar domain walls, it 
is suitable in particular for ferroelastic ferroelectrics such as crystals of the KDP family in 
which a dense pattern of 180° domains is known to exist. ll|, |l2j It was found that in a wide 
temperature range below the Curie point of crystals of RbH 2 P04 and deuterated KDP, the 
piezoelectric coefficient is greatly enhanced compared to its expected value for single 



domain samples. [I3j Recently, this was confirmed by simultaneous measurements of permit- 
tivity £3, elastic compliance and piezoelectric coefficient d^ of RbH^PC^. A thorough 
discussion of <i 36 in this case will be the subject of a forthcoming paper Kopal et al. Here 
we comment on the extrinsic part of £3. For simplicity, let us assume that = e c . This is 
not an unreasonable assumption since the assumed surface layer for KDP-type samples can 
be supposed to have a similar chemical composition as the bulk. Then the extrinsic part of 
Eq. (JSD reduces to 

B 21n2(l + £) 

e c R® q (e c + yje a e c ) 

Numerically, the second term in the brackets represents a small correction to the first term 
when approximations f^j and (HI) are valid. If it is neglected, we obtain as an approximation 



This shows that a very thin surface layer can lead to a considerable extrinsic enhancement of 



permittivity. Nevertheless, the simple implication: d —>■ =>- £ 2)C xtrmsic oo is not correct, 
because the assumptions needed for validity of (j3J), (HI), (J5J) and (jBJ) are violated if d is small 
enough. 

A more general formulation of the restoring force can be used to calculate the extrinsic 
part of d 3G for the same geometry of domains. In a recent paper [lj] Sidorkin deduced the 
dispersion law of wall contributions to permittivity, however, in his treatment the existence 
of a surface layer is not explicitly considered. 

It was shown beyond any doubt that small motions of 90° domain walls are responsible for 
a considerable enhancement of permittivity £3 and piezoelectric coefficients d$i, d^ in poled 
ferroelectric ceramics. Ua, ll6j One of the sources of the restoring force responsible for these 



wall contributions is the elastic energy at grain boundaries. [l7] Since these boundaries may 
differ in chemical composition from the bulk of grains lj| , surface layers can be expected to 
form so that the mechanism proposed in the present paper may also play a role in ceramic 
samples. 
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